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A model calculation is presented with the aim to study the interplay between magnetic and 
structural transitions. The model consists of an orbitally doubly degenerate conduction band and 
a periodic array of local moments. The band electrons interact with the local spins via the s-f 
interaction. The interaction of the band electrons with phonons is introduced by including band 
Jahn- Teller (J-T) interaction. The model Hamiltonian, including the above terms, is solved for 
the single particle Greens function. In doing this an ansatz for selfenergy of electrons, which was 
developed earlier has been utilized. The quasiparticle density of states (QDOS) and hence the 
orbital populations are calculated treating the ferromagnetism of local moments in the mean field 
approximation. The critical value of electron-phonon interaction (G) for the appearance of the band 
J-T distortion is higher in the ferromagnetic state. The strain appears at a critical temperature 
(Ts) when G is greater than the critical value. The onset of ferromagnetism at Tc {< Ts) arrests 
the growth of the strain. It is concluded that the magnetization hinders the structural transition. 
The quasiparticle density of states are presented to interpret these results. 

PACS numbers: 71.10.Fd, 75.30.MB, 75.30Vn 



I. INTRODUCTION 

Interplay of order parameters originating from differ- 
ent degrees of freedom of electrons is very relevant in con- 
densed matter. One of the recent interests in this area is 
related to the mutual influence of band Jahn- Teller (J-T) 
effect and ferromagnetic order. The former results from 
the removal of orbital degeneracy and the latter from the 
lifting of spin degeneracy. The magnetic order can origi- 
nate either from the same electrons that are involved in 
the band J-T effect or from the electrons in some other 
states. The first case corresponds to narrow band solids 
which are usually described by the two-fold degenerate 
Hubbard model and is applicable to intermetallic com- 
pounds. The interplay between the above two order pa- 
rameters had been considered earlier and it was shown 
that the magnetic order tends to suppress the sponta- 
neous distortion resulting from charge-lattice interaction. 
In the present work, we are considering a situation where 
the electrons involved in J-T splitting and the magnetic 
order are different. The doped manganites^'^ and Heusler 
alloys^'"*'^ are some of the real systems that represent this 
case. Doped manganites shows a rich variety of phenom- 
ena such as several forms of magnetic, orbital and charge 
ordering^'^ The interesting physics of manganites is due 
to the dynamics of the d-electrons of the Mn ion. In 
the host material LaMnO^, which is a Mott insulator, 
there are four electrons per AIn ion. The d-orbitals are 
split by the crystal field in to two, namely the triply 
degenerate (ignoring the spin) t2g states which are well 
localized and the doubly degenerate eg states which are 



extended and form a doubly degenerate band. Out of 
the four d-electrons, three occupy the t2g-states and the 
remaining electron is in the Cg-band. The Mn ion is in 
Mn'^+ state. Due to the large Hund's rule couphng in 
this system, the three <:2g-electrons have their spins ori- 
ented parallel to each other making a localized spin of 
S = 3/2. This is again strongly Hund's rule coupled to 
the Bg-electron. Since the Cg is an extended state, the 
electron in this state can hop from lattice site to lattice 
site. The hopping combined with Hund's rule coupling is 
responsible for the long range magnetic order that exists 
in this system. This is known in literature as the dou- 
ble exchange mechanism. Another way of looking at the 
system is that there is a localized spin at each lattice site 
and the band electrons interact with these spins via an 
intraatomic exchange interaction. This is known in lit- 
erature as the Kondo-lattice model (KLM). Similar sce- 
nario also prevails in Heusler alloy like Ni2MnGa, where 
Mn posses localized moment and the magnetic proper- 
ties of the host system can be understood from the KLM 
with the carrier concentration (e^-electrons) of one per 
Mn atom. When the trivalent rare-earth ion is partly 
replaced by divalent ions like Ca, Sr or Ba, the mate- 
rial exibits CMR properties undergoing transition from 
the paramagnetic-insulator state to the ferromagnetic- 
metal one. The localized spin of 5* 3/2 is retained but 
there is a decrease in the electron concentration in the 
eg-band. The carrier concentration is less than one per 
atom. While studying the CMR and the associated mag- 
netic and insulator-metal transitions, it was realized that 
KLM alone is not sufficient to understand the physics of 
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manganitcs. It is now accepted that eleetron-phonon in- 
teraction plays an important role^'^'^"'^^. These materi- 
als crystallize in perovskite structure and the Mn ion has 
Oq octahedron as its immediate environment and there- 
fore interacts strongly with the distortions in the octahe- 
dron. Obviously, this introduces J-T effect which lifts the 
degeneracy of the e^-band. The spontaneous distortion 
associated with the J-T effect exists when the lowering 
of the band energy is more than the increase in elastic 
energy due to the strain. The simplest way to describe 
the band J-T effect is to incorporate the interaction of 
the Cg-electrons with the lattice distortions (phonons) in 
the KLM model. When this is done, it is pertinent to 
examine the interplay between the band J-T effect and 
magnetism in the model system. There are experimental 
results supporting such an interplay in manganites and 
in Heusler alloys. It was observed that the J-T distor- 
tion in {La-Y)CaMn03^^ and {La,Pr,Ca)Mn03^^ is 
reduced in ferromagnetic state. The suppression of the J- 
T distortion in {Nd, Sm) 1/2 -Sr 1/2 AfnOs in ferromagnetic 
state under high magnetic field has also been reported^^. 
The coexistance of the J-T distorted phase and the fer- 
romagnetic phase have been reported in Rh2CoSn^ . J-T 
splitting exists in these systems below the transition tem- 
perature Tc- 

Therefore, in order to study the interplay of structural 
and magnetic transitions, we consider a model where the 
band electrons, which are approximated to be s-electrons 
in doubly degenerate extended states interact intraatom- 
ically with a periodic array of localized spins (a Kondo 
lattice). We provide for the spontaneous lifting of the 
degeneracy of the band states by including a band J-T 
interaction. The lifting of the degeneracy is signalled by 
the appearance of strain. The presence of the long range 
magnetic order is characterized by the nonzero value of 
the magnetization. We study, selfconsistently, the strain 
as a function of the J-T coupling constant for different 
carrier (cg-electron) concentrations and the dependence 
of the strain on temperature. In the latter case, the mag- 
netization of the local moments which is caused by the 
exchange interaction between the eg- and f2g-electrons 
(Kondo interaction) determines, decissively, the temper- 
ature dependence of the strain. 

II. MODEL HAMILTONIAN AND ITS 
APPROXIMATE SOLUTION 

The Cg-electrons moving in the doubly degenerate band 
are described by 
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Hs = {Tij — Il5ij) C^if^Cajo- 

a=l ij(T 

= Y.^ei'k) - f^)cli,^Cc,Ua (1) 

Tij is the hopping integral for hopping of the elec- 
trons from lattice site i to j. cl^ia-icaia-) is the cre- 



ation(annihilation) operator for an electron in the a-state 
on the lattice site i with spin a. a = 1, 2 is the band in- 
dex. ^ is the chemical potential. e(k) is the band energy 
related to by 

^^i = ^E^(k)e*-^''*-^^-' (2) 

k 

The band electrons interact with the localized spins 
via the intraatomic exchange interaction of the coupling 
strength J and this is described by 

Hsf = -J ^CTaJ ■ Sj 

= —-j''^^{Za-Sjnaja-+Sj C^aj-a^aja) (3) 

a is the spin of the band electron and S is the localized 
spin (total spin of the three t2g electrons). Uajcr is the 
number operator for the electron in the state a at the 
lattice site j with spin a. At the outset itself we assume 
a ferromagnetic interation (J > 0). is a sign factor, 
Za = Scri - 6a I and Sj = Sj + iZaS^. 

The electron density in the degenerate band couples 
to the static elastic strain through the J-T interaction. 
In the case of a tetragonal distortion, this interaction is 
described by^^.ie.i? 

HjT = Ge ^(nika - n2k<7) = Ge ^(nij^ - n2ia)- (4) 

k,(T itr 

G is the strength of the J-T coupling and e is the lattice 
strain given by 

G . 

^ = TVC^ Zl(< "'1''^ >-< n2ia >) (5) 

where Co is the elastic constant. It is clear that Hjt 
tries to create a difference in the occupation of the two 
degenerate bands. The difference in occupation leads to 
the building up of the strain. Thus, under suitable con- 
ditions, there is a spontaneous splitting of the bands and 
building up of strain which indicates a structural transi- 
tion. The building up of the strain however leads to an 
increase in the lattice elastic energy which is given by 

Hl = ^NCoe'' (6) 

Where A'^ is the total number of atoms. Since this term is 
a c-number and we are not looking for the ground state 
whose energy has to be minimum, we leave this term out 
of our consideration. Then the Hamiltonian of the model 
system we are considering is 

H = H, + Hsf+ Hjt. (7) 

The model Hamiltonian Eq(7) obviously cannot be solved 
exactly. However, in an earlier work^*, we have proposed, 
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for the Hamiltonian without the J-T term, an approxi- 
mation scheme, which is rehable in the hmit of low carrier 
concentration. We will exploit that scheme in solving the 
present model. Firstly, without resorting to any approxi- 
mation, we can absorb the J-T term into Hg by modifying 
the band energies for the two bands as 



e„(k) = e(k) + i-irCe. 



Then we have 



akcr 



and the total Hamiltonian is given by 
H = Hs + Hsf. 



(8) 
(9) 

(10) 



In order to calculate the strain ( caused by the structural 
transition), one has to calculate the one-electron Greens 
function 
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(11) 



That means, one has to calculate the selfenergy 
of the electron in the presence of Hgf. There are some 
exact results available for the selfenergy in certain lim- 
iting cases, namely, the zero band width limit for all 
temperatures^^ and the finite band width but ferromag- 
netic saturation (T = 0) limit^^'^^. In addition, using 
Mori formalism, the result for the second order pertur- 
bation theory is also available^^. We propose an ansatz 
for the self energy that reproduces the known limiting re- 
sults and in addition, satisfies the strong coupling limit. 
This can be taken care of by making a high energy ex- 
pansion and evaluating the first four spectral moments^^ . 
We thus have a selfenergy which fulfills a) the zero band 
width limit for all temperatures and coupling strengths, 
b) T = limit for all band widths and coupling con- 
stants, c) the weak coupling limit (OJ^) and d)has the 
correct high energy behaviour. The ansatz for T,aa{E) is 
given by 

^aa[E) - --Jm„ + -J — T- r-. (12) 

2 4 1 - b^G ao [E - ^Jnna) 



Here 



g..(e)=ix:g,.(e)= 1^:^-1^ (13) 



rucr = Za (S^). The ansatz assumes a k— independent 
selfenergy. As Hgf is a local interaction, the energy de- 
pendence of the selfenergy is the deciding factor in rela- 
tion to the electron density of states. The parameters a„ 
and bcr and are fixed by rigorous high energy expansions 
to fulfill the first four spectral moments: 
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J (14) 



It should be mentioned that this ansatz is valid only in 
the limit of low carrier density. Since we are interseted in 
simulating systems with low carrier density, it is justified 
to use the above ansatz. From GaWaiE) one can obtain 
the spectral density Sa]ia{E) and the density of states 
Paa{E) from the well known relations 



Saka{E) = --ImGaka{E) 

PaaiE) = 1^5„k.(i;) 



(15) 
(16) 



From the knowledge of the density of states, the expec- 
tation values can be evaluated: 



(n„,) = / dEf_{E)p„„{E) 



(17) 



Where f^{E) = 1/ (l e'^^) is the Fermi function with 
/? = 1/kT. The chemical potential is fixed by the 
constraint 



naa) = constant. 



(18) 



For a given set of the model parameters n, J, S and 
G and for a fixed T, the occupencies of both lower and 
upper sub-bands for each spin directions is computed sclf- 
consistently. After having the self consistent solution, the 
lattice strain 



G 



(19) 



is calculated. The average occupation of eg-orbitals (naa) 
can be numerically obtained using a model density of 



states for the "free" Cg-band: 



po{e) = aJi 



E 








In 


E^ 


D 



(20) 



Where ^4 is a normalization constant and D is half the 
width of free Bloch band. In order to calculate {naa), we 
require {S^), since this enters into the selfenergy. The 
local moment system is described within the mean field 
approximation and is represented by Brillouin function. 
The effective field seen by the local moment is deter- 
mined by mutual exchange interaction which fixes Tc- 
Numerical results are given where all the energy parame- 
ters are normalized in terms of the free bandwidth {2D). 
Ideally, one should get the magnetization (S^) selfcon- 
sistently out of the calculation. However, it is a very 
involved problem. Therefore, we treat it as a parameter 
and obtain its value at any temperature from the Bril- 
louin function assuming a value for Tc- 

Using the above theory, the results for the density of 
states and the strain are presented in the next section. 
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III. RESULTS AND DISCUSSION 

First we consider the T = case for two extreme 
situations of the local magnetization, namely, the para- 
magnetic {{S^) = 0) and the saturated ferromagnetic 
{{S^) = S) state. In Fig.l we have plotted the strain as 
function of the J-T coupling constant G. We find that 
unless the value of G exceeds a critical value, there is no 
spontaneous splitting of the orbitally degenerate bands. 
The splitting of the bands takes place and the lower of 
the split bands is occupied more by the electrons than 
the upper one in order to lower energy. This is energet- 
ically favoured only for a sufficiently large G. When G 
is further increased, the lower of the split bands is more 
populated and therefore the strain increases as shown in 
the figure. Though not shown in the figure, the criti- 
cal value of G also depends on the carrier concentration. 
Larger is n, smaller is the critical value of G. Another 
feature which is displayed in Fig.l, is the role of the mag- 



0.4 




FIG. 1: Lattice strain as a function of electron-phonon cou- 
pling constatnt G at T = with and without localized mag- 
netization, n = 0.7, J = 1, S = 3/2 and Co = 1. 

netization (S^). Critical value of G to induce J-T effect is 
much higher for the saturated local magnetization situa- 
tion as compared to that of the paramagnetic situation. 
With increase of G, it is observed that the gap between 
the two curves decreases. The difference in the critical 
value of G appears to be a result of the competition be- 
tween the mechanisms leading to the lifting of orbital 
and spin degeneracies of the e^-state. In the ferromag- 
netic situation, the effective field lifts the spin degeneracy 



so that stronger electron lattice interaction is necessary 
for the strain to appear. As the strain becomes larger, 
the energy gain by the J-T distortion dominates over the 
energy gain by lifting of the spin degeneracy. So the in- 
fluence of magnetization becomes less important. The 
interplay becomes more interesting when the energy gain 
due to the two mechanisms is comparable. 

In order to examine this interplay in more detail, we 
study the temperature dependence of the strain in Fig. 
2. The structural transition temperature is the tem- 
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FIG. 2; Lattice strain as a function of temperature for various 
values of Tc (curves a, b, c, d and e). n — 0.7, J = 1, S = 3/2, 
Co ~ 1 and G = 0.42. Brillouin function is also plotted (curves 
1,2,3 and 4) in the same figure as a guide for the eye. 

perature at which the strain goes to zero. We choose 
the parameters such that Tg = 800K. Then we study 
the effect of magnetization on strain by varying Tc such 
that i) Tc < Ts, ii) Tc « Ts and iii) Tc > Ts. The 
Fig. 2 displays the T-dpcndcncc of the strain(curvcs a, 
b c ,d and e) and the magnetization (which is a Brillouin 
function for a given Tc) for different Tc's (curves 1, 2, 3, 
and 4). Whatever is Tc, at T = 0, (S^) = S. Therefore, 
the effect of the magnetization on the strain is indepen- 
dent of Tc and leads to a maximum decrease of strain. 
As Tc increases, the strain goes up with increasing T 
and the rate of increase is higher for lower Tc so long 
as Tc < Ts- therefore, there always appears a peak in 
the curves. When Tc « Tg, the peak is very faint and 
the strain becomes very small well before Tg is reached. 
For Tc > Ts, the strain is nonzero only when T is much 
lower than Tc (curve e of Fig. 2). 

It is clear that at T = 0, the presence of magnetization 
causes a redistribution of electrons between the orbital 
levels by creating a population difference between the 
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spin levels. Such redistribution is the cause of suppres- 
sion of strain. As T increases, the spin level occupancies 
tend to equalize and that polarises the orbital levels fur- 
ther resulting in the increase in the strain compared to its 
T = value. Since for smaller Tq the magnetization de- 
creases faster with increasing T, the increase in strain is 
also faster. For Tc > Tg, due to the choice of parameters, 
as expected, there is no strain between Ts < T < Tc- 
When T is much lower than Tc, the occupancy of the 
spin levels is stabilized and the system can lower energy 
by further redistribution of electrons between the orbital 
states. That is why the strain is finite for T much less 
than Tc (curve e of Fig. 2). The results obtained are 




Energy Energy 



FIG. 3: Quasiparticle density of states as a function of enr- 
egy for various values of localized magnetization (S^) in the 
absence of Jahn- Teller distortion i.e for G = 0. Full line for 
spin up and dotted line for spin down. J = 1 and S = 3/2 

interpreted with the help of the quasiparticle density of 
states (QDOS). Before discussing the results of the full 
problem, to fix up a reference for further discussion, we 
want to present in Fig. 3, the QDOS for the case of a pure 
KLM, that is, for the case of G = 0. The QDOS consists 
of two subbands for each spin direction separated by an 
energy of the order of ^J{2S + 1). The separation of the 
bands is independent of T but the spectral weights of 
these subbands, however, depend on T through the value 
of (5'^). For example, at T = ((S'^) = S), the spectral 
weight of the upper subband for t-states is zero. The rea- 
son for this is easy to understand. At T = 0, the local 
moment system is saturated. Therefore, for an |-electron 
there is no chance to spin-flip by involving a correspond- 
ing spin- flip of the local moment system. That means, at 
T = 0, as far as the t-electron is concerned, only the Ising 
part of Hgf operates resulting simply in a rigid shift of 
QDOS. As we see from Fig. 3, the spectral weight of the 
J,-states in the lower sub band is finite. This is because, 
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FIG. 4: Quasiparticle density of states (for lower sub- 

band(Q = 1) in the positive half of the frame and for upper 
sub-band(a = 2) in the negative half of the frame) as a func- 
tion of enregy for various values of electron-phonon coupling 
constatnt G. Full line for spin up and dotted line for spin 
down. Thin vertical line shows the position of the chemical 
potential, n = 0.7, J = 1, S = 3/2, Co = 1, (S"^> = S and 
T = 0. 



for a J.-electron, even at T = 0, spin-flip is possible. Fur- 
ther more, when a J,-electron flips its spin, it lands as 

an t-electron. Therefore the nonzero QDOS of the |- 
electron should be in the same energy region as that of 
t-electron. For a t-electron there is another possibility. 
It can have repeated magnon emission and absorption. 
That is, in a sense, it propagates in the lattice dressed 
by a cloud of magnaons. This is a stable quasiparticle, 
which we call the magnetic polaron. Obviously, at T = 0, 
there is no possibility of magnetic polaron for t-electron. 
As T increases {{S^) decreases), the spin flip processes 
are allowed for both the spin directions and therefore the 
spectral weights in both the subbands are nonzero for 
both the spin directions. At T = Tc ((S"') = 0), the 
spectral weights of t- and t-states in the two subbands 
become equal as it should be. We note the assymetry 
with respect to the centre of the free band. This orig- 
inates from the renormalization of the atomic levels by 
the s-f interaction^^. 

Now we consider the further splitting of these bands 
due to the J-T effect. That is, when the degeneracy of 
the Cg band is lifted due to the J-T effect, each of the 
subbands of Fig. 3 for each spin direction discussed above 
is again split into two as shown in Fig. 4. Noting the po- 
sition of the chemical potential, we see that the upper 
subbands for both values of a can be ignored since they 
are never populated. Only we have to keep in mind the 
change in their spectral weights. From now on, when we 
speak of subbands they are the two eg bands, correspond- 
ing to a = 1, 2 as displayed in the positive and negative 
halves of the frame. 

We want to understand the dependence of the strain 
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on the coupling constant G and the influence of (S^). In 
Fig. 4, we take the case of saturation ((5^) = S). When 
G is less than the critical value, the two Cg bands are 
degenerate. As G approaches the critical value of 0.42, 
there is a slight splitting of the two subbands and the 
diflference in their population becomes non zero. At the 
same time, it should be noticed that the position of both 
the subbands shifts to lower energy, so that, on the whole, 
the energy of the system is lowered by the splitting. For 
a small increase in G, the splitting increases and the in- 
crease in strain is very large. Further, the shifting of the 
subbands to lower energy is also large so that the energy 
of the system is much lower. Any further increase in G 
does not have much effect on the quasiparticlc spectrum 
any more which means the strain saturates. It should be 
emphasized that the strain is not introduced by hand but 
comes selfconsistently out of the model. 
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FIG. 5: Quasiparticle density of states(for lower sub-band in 
the positive half of the frame and for upper sub-band in the 
negative half of the frame) as a function of enregy for various 
values of localized magnetization (S^). Full line for spin up 
and dotted line for spin down. Thin vertical line shows the 
position of the chemical potential, n = 0.7, J = 1, S = 3/2, 
Co = 1, T = and G = 0.42. 



In Fig. 5 we now fix the value of G at 0.42, so that the 
model prefers the J-T split situation and look at the in- 
fluence of (5^). At saturation {{S^) = S), the spectral 
weight of "f-subband is such that the a = I subband is 
slightly more populated compared to the a — 2 subband. 
As a result, a small strain appears. As {S^) decreases, 
this spectral weight is modified in such a way that the 
occupation of the a = 1 subband is more than that of the 
a = 2. This results in an increase of the strain with T as 
depicted in Fig. 2. This trend continues until (S^) = 
and at that point, the strain is as if there is no exchange 
interaction in the model. In short, the effect of magneti- 
zation is to hinder J-T splitting. 

In order to examine the interplay of magnetization 
and the strain as a function of T, it is necessary to in- 
clude the temperature variation of both the strain and 
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FIG. 6: Quasiparticle density of states(for lower sub-band in 
the positive half of the frame and for upper sub-band in the 
negative half of the frame) as a function of enregy for various 
values of temperature T. Full line for spin up and dotted line 
for spin down. Thin vertical line shows the position of the 
chemical potential, n = 0.7, J = 1, S = 3/2, Co = 1, and 
G = 0.42 . 



the magnetization. Therefore, we consider a specific case 
of Tc = 250K and = 800K and the corresponding 
QDOS at different temperature are displayed in Fig. 6. 
Starting from T — 0, the a = 1 subband occupation in- 
creases as T increases upto Tc- Above Tc, the situation 
is reversed and both the Q-subbands are equally occupied 
at Tg. Therefore, the maximum of strain occurs at Tc as 
shown in curve (b) of Fig. 2. 



IV. CONCLUSIONS 

The model Hamiltonian is solved by first absorbing the 
J-T term into the band electron term and then utiliz- 
ing an interpolation ansatz for the sclfencrgy. The band 
splitting and through it the strain in the lattice due to J- 
T interaction has been determined selfconsistently. The 
strain as a function of the coupling constant at T = is 
studied with and without the presence of local moment 
ordering. It is found that a minimum value of G is re- 
quired for the strain to appear. This critical G is larger 
if the local moments are ordered or alternately if there is 
an external magnetic field. The temperature dependence 
of the strain is studied by assuming different values for 
the magnetic transition temperature. The study indi- 
cates that there is a strong interplay between the mag- 
netic and structural transitions. It is observed that the 
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growth of the strain appearing at T, (> Tc) is arrested 
with the onset of ferromagnetism and tends to a lower 
value determined by the magnetization at T = 0. This 
means the removal of spin degeneracy is not conducive 
to the removal of orbital degeneracy. The results are ex- 
plained on the basis of the QDOS. The basic ingredients 
of the model are the band J-T effect, long range mag- 
netic order and their mutual interaction. The possibility 
of different hopping between the degenerate bands is not 
considered for simplicity and therefore the results cor- 
respond to large J-T effect. The inter-orbital hopping, 
which is sometimes considered, would remove partially 
the degeneracy of the state. Also the magnetism which 
should evolve within the model is treated as a parameter. 
The detailed comparison with experiment is therefore not 
attempted. However, the general trends of the results 
related to suppression of J-T strain and coexistance of 
two phases are in tune with the experimental^^'^^'^^ ob- 
servations in a Ca-doped manganite and Heusler alloy. 
The calculation is based on the ansatz used for sclfen- 
ergy which is valid only in the limit of low charge carrier 
concentration. If this needs to be relaxed, the interaction 
among the band electrons has to be taken into account 
and naturally the ansatz for the selfenergy has to be mod- 
ified. 
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